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Abstract 

We generalize the geodesic rule to the case of formation of higher codimensional 
global defects. Relying on energetic arguments, we argue that, for such defects, the ge- 
ometric structures of interest are the totally geodesic submanifolds. On the other hand, 
stochastic arguments lead to a diffusion equation approach, from which the geodesic 
rule is deduced. It turns out that the most appropriate geometric structure that one 
should consider is the convex hull of the values of the order parameter on the causal 
volumes whose collision gives rise to the defect. We explain why these two approaches 
lead to similar results when calculating the density of global defects by using a theorem 
of Cheeger and Gromoll. We present a computation of the probability of formation of 
strings/vortices in the case of a system, such as nematic liquid crystals, whose vacuum 
is RP 2 . 
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1. Introduction 



Topological defects [1] are static solutions of the Euler-Lagrange equations of a 
classical field theory, whose non-triviality is guaranteed by the homotopy type of the 
vacuum. The stability of a topological defect can be ascribed to the infinite amount of 
energy required to reduce it to the trivial configuration, whose energy is zero [1]. This 
highlights one of the reasons for the methodological importance of topological defects 
in field theories: they provide a glimpse of non-perturbative aspects of the vacuum of 
a model, a vacuum for which very few things are generically known. 

The geodesic rule [2]-[ll] is a way that allows us, in principle, to make predictions 
about the density of such topological defects. There seems to be a significant body 
of experimental [12]- [16] and numerical results [17], [18] that makes us be reasonably 
confident about the validity of the geodesic rule for the case of global defects. The 
formulation and the validity of this rule for local defects, however, is a much more 
contentious issue, both experimentally and numerically [19], [20]. For this reason so we 
will confine our analysis to just global defects in the present work. It is worth noticing 
that the geodesic rule is not the only known way for producing topological defects. 
Another way is through order parameter oscillations [21], [22], [23] whose contributions 
are dominant is when there is no significant dissipation during bubble collisions. 

To be more concrete, we will assume that the phase transition after which a topo- 
logical defect is formed is of first order and proceeds by bubble nucleation [1], [10] ,[11]. 
Then the geodesic rule states that, when two such bubbles collide, the order parameter 
(spacetime scalar) <fi smoothly interpolates between its values in the bubbles 0i, 02 
following the shortest path (geodesic) on the vacuum manifold. There are two ways to 
justify the validity of geodesic rule: The first relies on minimizing the energy density 
of the model under consideration in the spontaneously broken phase [4], [5]. The sec- 
ond relies on the stochastic nature of the bubble collisions in the broken phase and the 
Markovian character of the bubble effective dynamics during their expansions and their 
collisions [24], [25]. In this paper we extend the arguments of both of these approaches 
to the case of formation of higher co-dimensional defects (monopoles, textures etc). We 
justify the reason why these two approaches lead to the same estimate for the density 
of global defects by using a theorem of Cheeger and Gromoll [26]. We use an integral 
geometric argument to determine the probability of formation of global strings/vortices 
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formed in a system whose vacuum is MP 2 as was also done in an earlier work. 

2. The energetic approach and totally geodesic submanifolds 

A classical, quantum or thermal field theory describing the phase transition that 
gives rise to topological defects, has as basic kinematic variables the order parameter(s) 
0, which, for simplicity, we will assume to be spacetime scalar fields in this work. Let 
the spacetime, on which the model is defined, be M fc , endowed with the usual flat 
metric, for concreteness. The scalar fields are sections of a vector bundle with base 
M. k and typical fiber an, irreducuble most of the times, representation of a compact, 
semisimple Lie group G, which incorporates the internal symmetries of the system. 
The vacuum then consists of all such maps that minimize an appropriate energy 
functional. In most cases we choose to focus our attention on the image of the space- 
time points under 0, rather than on the set of such maps themselves, as the basic 
kinematic variables. Such a choice is largely unimportant however, especially within the 
domain of perturbation theory around a smooth 0, like the ones used in the Ginzburg- 
Landau description of phase transitions, which we will be working with in what follows. 

We assume that the vacuum of the model is an m-dimensional Riemannian manifold 
A4, endowed with a metric g. In most cases of interest, as was mentioned above, 
M. is a semisimple, compact Lie group G, endowed with its Cartan-Killing metric 
g = Tr(u~ l du) 2 , u e G. Then the Cartan-Killing metric is non-degenerate and 
positive definite, so (G, g) is, indeed, a Riemannian manifold. The vacuum M. is 
also very frequently a homogeneous space G/H, endowed with the induced metric 
from g. The energy functional of the model can be written as 

m = [ da0x^)+y(0)} (i) 

where a — 1, 2, . . . , k and V(4>) denotes the potential energy density of the model. 
The possible existence of a mass term ^0 2 is assumed to be incorporated inside 
V(4>). Local energy extremization, under an infinitesimal energy variation 5E gives 

8E[<t>] = [ {(d^ida^ + SVm (2) 

All the points of the stable vacuum have the same potential energy density, so we can 
set 5V(4>) =0 on M. Then the vanishing of the first variation of the energy reduces 
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to 

dj<j) = (3) 

This implies that the energy between two points of the vacuum <pi and 02 wm be 
extremized when <50 is constant [4], [5]. This constant can be made equal to zero by 
an affine reparametrization. Such a result can be heuristically interpreted, as stating 
that the scalar field must follow a path of "stationary length", with respect to per- 
turbations, when it is arc-length parametrized, namely it must follow a geodesic, with 
endpoints 0i and 02. Such a geodesic is locally length minimizing [26], [27], as is 
also noted right after eq. (28) below. Generally, the variation 5 is in the space of all 
maps that are sufficiently smooth for our purposes. However, since we have decided 
to focus on the image of points of M fc under 0, rather than on themselves, the 
variations 5 effectively reduce to variations on M. 

One can extend, still heuristically, this argument to the case of higher co-dimensional 
defects. For concreteness, let us consider the case of monopoles, which are pro- 
duced by the collision of four causally disconnected bubbles with scalar field values 
0i, i — 1, ... ,4. According to the geodesic rule, will interpolate between each pair 
of bubbles, along the geodesic(s) 0^-, i, j = 1, . . . , 4, i ^ j that join the values of in 
these bubbles, on M. Consider a 2-simplex in M with vertices 0j and edges 0^-. 
By analogy with the geodesic rule, and as we will prove in what follows, we should 
require the faces of this simplex to be least-area manifolds having as boundary the 
union of the geodesic segments (fiij. We can then determine the homotopy properties 
of the 2-simplex, and decide, in particular, whether it corresponds to a trivial element 
of the second homotopy group 7r 2 (.M). If it does, then a monopole will not be formed, 
on topological grounds, from this bubble collision. Otherwise, a monopole formation is 
possible. By enumerating the number of simplices with trivial and non-trivial images 
in n 2 (M.) with vertices uniformly distributed on A4, we can provide an estimate 
of the number of monopoles that will be formed in this phase transition. One can 
inductuvely extend this construction to the case of higher co-dimensional defects. 

We want to derive a local analytic condition that, if possible, is satisfied by the 
heuristic multi-dimensional extension of the geodesic rule described in the previous 
paragraph. Such an extension applies to the formation of n-codimensional defects. The 
homotopy group guaranteeing the stability of such defects, when non-trivial, is 7ik- n -i- 
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As an example of this homotopic classification, when k — 3 and n = we get the 
condition for the stability of monopoles expressed by the non-trivial n 2 . Similarly, 
we get conditions for the stability of strings/vortices, textures and other defects. Here 
n is the dimension of an appropriate sub-manifold Af properly embedded in M. 
as I : Af — > M.. Although / is a proper embedding, on physical grounds, i.e. no 
self-intersections of Af should be allowed, it turns out that such a condition is not too 
overly restrictive. Indeed, the following computations can be appropriately generalized 
when I are immersions with finite self-intersections, resulting in integral fc-varifolds 
for instance, without any major geometric modifications [27]. We allow for Af to 
have a piecewise smooth boundary dAf. This boundary is composed of simplices 
iteratively generated by the values of fa, i — 1, . . . , n + 2 in the colliding bubbles, 
the geodesies faj, i, j — 1, . . . , n + 2, i ^ j joining them, and their higher dimensional 
counterparts. Let g = l*g indicate the induced metric on Af and x G Af. We 
choose an orthonormal basis e\, . . . , e n , v n +i, . . . , v m of T X A4, with e^, i — l,...,n 
ortho normal in T x Af and — n + 1, . . . , m othonormal in the fiber N x Af. Then 
any X G T X M can be decomposed as X = X T + X 1 - where X T e T x Af and 
X 1 - E N x Af. Let V and V denote the Levi-Civita connections on M and Af, 
respectively. Then, for X, Y e TAf [28], 

VxY = (V X Y) T (4) 

Consider a smooth isotopy ("smooth deformation") of Af which keeps dAf fixed. 
Such an isotopy is a map 

F(x, t) : Af x (-e, e) -> M (5) 

with e > 0, which a smooth enough diffeomorphism F(-,t) : Af — > Af, for each 
t G (— e, e). The fact that such a diffeomorphism exists is guaranteed by the implicit 
function theorem. The vector field generating F(x,t) is 



dF(x,t) 



dt 



(6) 



t=o 



and let the induced metric on F(Af,t) be denoted by g t . Obviously g = g. 
According the area formula [29] , 

Vol{F(Af,t))= [ J(x,t)d n x (7) 
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where J(x,t) denotes the Jacobian associated with the action of F t . This Jacobian 
can also be expressed in terms of the metric as [30] 



J(x, t) = a/ det 



9t 



which means that 



d 



-Vol(F(Af,t)) 



t=o 



d / det gt 
'tf\dt\ det g 



a/ det g d n x 



(8) 



(9) 



Since g t is a positive-definite symmetric matrix, it satisfies the identity 

detg t = exp(tr\ng t ) 

where tr denotes the trace of g t and we get 
d 



dt 



^/det 



9t 



t=o 



(10) 



(11) 



which implies, taking into account that {e,}, i = l,...,n is an orthonormal basis, 
that 

d 



det g t 
dt V det go 



t=o 



1 . dg t 
- tr — — 

2 dt 



(12) 



t=o 



where we have used that the operations of tr and differentiation of g t commute with 
each other. Because V is a Levi-Civita connection, we have 

d 



dt 



g t {e h ei) = 2g(y Ft e h ei) 



and because it is torsion- free, namely 

V Ft ei- V ei F t = [F t , ei ] =0 

we find 



d 

dt 9t{ei ^ ] 



£(V ei F t , ei ) 



(13) 

(14) 
(15) 



t=o 



which, upon using Leibniz's rule with the Levi-Civita condition once more, gives 

V e Xg{F u e i )\ = ~g{V ei F t ,e i ) + ~g{F t ,V ei e i ) (16) 

resulting in 

d I det g t r _ /r , _ ^ ^ ^ ^ 



dt y det go 



= J2VeA~9(Ft,e l )}-~g(F t ,V ei e l ) 
t=o i=1 
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Because g(F t , ei) is a scalar field and by the definition of divergence [27], [30] 

n 

div Ar X = Y,9(Ve i X,e i ), VlGT,Af (18) 



i=l 

we find that 



d / det g t 



= div N F? -J2g(Ft,^ ei ei) (19) 
t=o i=1 



(it y det 

Combining Stokes' theorem and the initial assumption that F t = on dN , we 
conclude that 

/ divrfF? y/detg~ d n x = (20) 

The second term on the right hand side of (19) can be re-expressed, due to linearity, 

as 

n n n 

~9(Ft, V^) = 9(F t , (V ei e i ) ± ) + £ g(F t , (V e ^) T ) (21) 

i=l i=l i=l 

We observe that the second term reduces to 

n n 

~9(Ft, (V ei e,) T ) = 9(Ft, V ei eO (22) 
i=i i=i 

By choosing a normal coordinate system around x G A/", we see that V ei ej = and, 
since this is a vector relation, it is true in all coordinate systems. Therefore 



J2~9(Ft,(Ve i e l ) T )=0 (23) 
i=i 

The second fundamental form at x E Af is a vector in N X M given by 

S(X,Y) = (V X Y) ± , VX,YeT x N (24) 

We have used the sign convention of [28] in this definition. The trace of the second 
fundamental form is called the mean curvature vector H, so, 

1 " 



H=-J2S(e l ,e t ) (25) 

i=l 

Then, the first term of (21) becomes 

n 

Y,~g(F t ,(V et e t ) ± )=n~g(F t ,H) (26) 



i=l 
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so, (9) eventually gives 



d_ 

di 



Vol(F(Af,t)) = I n g(F t , H) y/detg d' 1 



t=o 



m 



X 



(27) 



Since the isotopy F(x, t) is arbitrary, we conclude from (27) that the multidimensional 
generalization of the geodesic rule that we are seeking, amounts to determining all N 
for which 



Sub manifolds for which H = are called minimal [31], [32]. Obviously, geodesies are 
minimal 1-dimensional submanifolds of A4. The word "minimal" may appear to be a 
misnomer, since all that we have determined are the extrema of the volume functional, 
which are not necessarily minima. This is indeed true globally. To determine the 
volume minimizing minimal submanifolds, one would have to use the second variation 
formula [28], [29]. Locally, however, a minimal manifold without singularities is volume 
minimizing [27]. To be more concrete, assume that x is not on the boundary dAf. 
It turns out, then, that there is a small enough open set U C A~, with x £ U such 
that 



when Af is minimal. It is worth pointing out that although an appropriately modified 
version of (27) is true for manifolds with singularities, as mentioned above, (29) does 
not hold in such a case [27]. An obvious class of embeddings / satisfying (27) are the 
ones for which S = 0, instead of the weaker H = 0. Such submanifolds Af are 
called totally geodesic [28] , because they have the important property that a geodesic 
7(5) of A4 starting at x E Af with an initial direction g T x Af always stays in 
Af. The order parameter of the broken phase in the colliding bubbles and the minimal 
geodesies joining them should be entirely inside Af. In any other case, there is at 
least a part of a geodesic segment which is in the complement of A", thus failing to 
extremize the energy (1). Such cases are ruled out by having requiring the extremiza- 
tion (2). This argument demonstrates that we must choose as physically relevant only 
the subset of minimal submanifolds that are totally geodesic. These are the higher 
dimensional generalizations of the geodesic rule that we were seeking. Obviously the 
geodesies themselves fulfill this requirement. 



H = 0, VxeAl 



(28) 



Voi(Af nu)< Voi(F(Af n u, t)) 



(29) 



3. Stochasticity, convexity and equivalence of results 
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In previous work [24], [25] one of us presented an argument leading to the geodesic 
rule on stochastic grounds. The argument relied, in an essential way, on the existence 
of two widely distinct time scales, one relatively short for the bubble collisions r c , and 
a much longer one for the bubble coalescence tm- With the additional assumptions 
of the Markovian character of the bubble dynamics and the almost isotropy of the vac- 
uum A4, a diffusion equation governing the evolution of the order parameter inside 
the bubbles in the low-temperature phase was derived. The short "time" asymptotics 
of the corresponding solution gave rise to the geodesic rule [24], [25]. 

There is no obvious a priori way to generalize this result and determine the higher 
dimensional analogues of the geodesies, call them Q, in the stochastic approach. 
However we can proceed inductively as follows: Let's consider the value of the order 
parameter fa, * = 1? 2, 3 inside three bubbles that collide to potentially give rise to 
a string-like defect. Assume that these three values are within the injectivity radii of 
each other so that the geodesic rule as argued in [24] holds. Let faj, i,j = 1, 2, 3, i ^ j 
denote the minimal geodesies joining the bubbles whose order parameter is fa and 
4>j. Due to the proximity of fa and fa there is just one such segment (mini- 
mal geodesic) for each pair ij, i ^ j. Let p G faj and q G (fiik with i ^ k. 
Then, the unique geodesic fa q must be contained in the set Q, due to continuity. 
Continuity also demands all pairs of interior points of Q be connected by segments. 
All the points of such segments must belong to Q. In Euclidean space we would 
call a set Q satisfying such a requirement "convex". We must be careful, however, 
since there are several inequivalent concepts of convexity in the Riemannian context, 
all of which coincide in the Euclidean case [23], [25], [27]. To state the problem: what 
we are seeking is a set Q C M, such that for any two p,q G Q there exists a 
segment fa q C Q such that fa q is the unique segment connecting p and q in 
M. The set Q satisfying this condition is called "convex" even in the Riemannian 
setting [26], [28], [30]. The geometric structure Q we wanted to determine, therefore, 
is the convex hull of fa, i = 1, 2, 3 in M. We see that this result can be inductively 
generalized for all higher dimensional topological defects, without any further difficulty 

The class of classical vacua M usually employed satisfies one further condition. 
As was pointed out above, M. is very frequently a Lie group (G,g), where g is 
the Cartan-Killing metric or one of its homogeneous spaces G/H with the induced 
metric. For such cases more can be stated. Indeed, let W, Z,V G TG be left- 
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invariant vector fields on G and [W, Z] be their commutator. If \W A Z\ = 
g(W, W)g(Z, Z) - [g(W, Z)f and the Riemann tensor is 

R(W, Z)V = V W V Z V - V Z V W V - V [W:Z] V (30) 

it turns out that the sectional curvature 

g(R(W,Z)Z,W) 



can be computed to be 



K{wz) = m,w\) (31) 

v ' 4\W A Z\ v ; 



Therefore any Lie group (G, g) has non-negative sectional curvature. The same is true 
for its homogeneous spaces G/H with their induced metrics. To see that, it suffices 
to notice that the principal fibration H — > G A G/H, with ix being the natural 
projection, is a Riemannian submersion. Then, a formula [33] gives 

KfaW, n*Z) = K(W, Z) + ^ g([W, Z^, [W, Z^) (32) 

where W, Z are in addition assumed to be orthonormal and, in accordance with 
the notation of Section 2, W 1 - denotes the normal component (fiber tangential) of 
W E TG. Therefore we are interested in the convexity properties of subsets Q of 
vacua M. which are compact manifolds of non-negative sectional curvature. Conse- 
quently the Ricci curvature along a particular direction of TA4, which is the average 
of the sectional curvatures along all two-planes containing this direction [28], [30], is 
strictly positive. The most important point for the stochastic approach isn't so much 
that the Ricci curvature on M. is positive, but that it has a lower bound. Such a lower 
bound guarantees that the diffusion equation that is used to establish the geodesic rule 
on M. has several "reasonable", from a physical viewpoint, properties such as the 
total conservation of heat (stoachstic completeness), uniqueness of the associated heat 
kernel etc [34]. Such properties can be expected to be true on general physical grounds 
even if the underlying smooth structure of M. were to be abandoned, a fact toward 
which the quantum/thermal cases seem to be pointing, as will be explained in the next 
Section. 

It is of some interest to generically ascertain that such convex Q, as predicted by 
the geodesic rule, exist for physically relevant vacua M. If this were not true, then 
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topological defects would not form even when they should have, despite the homotopic 
arguments and the experimental and numerical evidence to the contrary. Then the 
stochastic approach, as we have developed it, would be invalid. As a crude first attempt, 
one can at least try to guarantee the existence of convex subsets of M. by estimating 
their maximum "size". An upper bound for the linear dimensions of such a convex 
subset is provided by the convexity radius [30] Conv(.M). It is defined for x G M. 
by 

Conyj, = sup {p : B x (r) is convex for all r < p} (33) 

where B x (r) stands for the solid ball centered at x of radius r. Then it can be 
proved that [35] 

Conv(^)<^M (34) 

where mj(A4) indicates the injectivity radius, which is the minimum distance around 
any x G M. for which the exponential map is a diffeomorphism. There are various up- 
per and lower bounds for the injectivity radius of a manifold, depending on conditions 
that the curvature, volume and other geometric characteristics of the manifold satisfy, 
bounds that we omit since are not explicitly needed in the sequel [36]. It is encouraging 
that the convexity radius is not identically equal to zero, although this may happen 
for some x G M. It would be useful to be able to find a lower bound estimate for 
the convexity radius. Such an estimate, however, is more difficult to obtain [36], but 
careful examination shows that this is not as serious a problem as it might appear at a 
first glance. Indeed, the existence of convex Q is guaranteed and explicitly provided 
by the interior construction of Cheeger and Gromoll [26] that leads to the proof of their 
structure theorem mentioned in the next paragraph. 

As we have seen, the generalizations of the geodesic rule resulting from the energetic 
and the stochastic approaches are embedded totally geodesic submanifolds and convex 
subsets of M, respectively. The question that naturally arises is whether these two 
classes of objects give rise to the same prediction about the density of topological de- 
fects. We expect that they should, since they describe the same physical phenomenon 
from two different viewpoints. This expectation turns out to be correct, as proved by 
Cheeger and Gromoll in the fundamental [26]. In that work, it was demonstrated that 
if Q is a closed connected convex subset of a Riemannian manifold Ai, then Q has 
the structure of an embedded n-dimensional submanifold of M. with a smooth totally 
geodesic interior and possibly non-smooth boundary dQ. Comparing the properties 
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of M in Section 2 and of Q of the present Section, we see that the Cheeger-Gromoll 
theorem implies that for each Q there is exactly one admissible M . Therefore, the 
predictions about the density of topological defects derived by the energetic and by 
the stochastic arguments coincide, as they should. The proof of the Cheeger-Gromoll 
theorem is attained in several steps and it is quite geometric. The methods employed 
as well as the constructions leading to the result, although highly transparent, appear 
to have only limited utility for physical purposes, so we skip them and refer to the 
original [26] for further details. 

4. Further comments and a sample calculation 

In this Section we are making three mutually loosely connected comments, related 
to the topics discussed above. First, we comment on the time scale in which the 
geodesic rule was established. As mentioned in Section 3, the stochastic derivation 
of the geodesic rule relied on the assumption that t m 3> r c [24], [25]. In the case 
t m ~ T Ci the strong mixing of the order parameters of the colliding bubbles suppresses 
the number of topological defects that can be formed. This can occur, for instance, 
when there is a high rate of bubble nucleation. To enforce this suppression explicitly, 
we can introduce as a regulator a "mass" parameter M in the diffusion equation that 
established the geodesic rule. Then the diffusion equation would become 

^ = (D + M 2 )V 2 (35) 

where D is an effective diffusion constant characterizing the dynamics of in the 
broken symmetry phase carried by the colliding bubbles. The exact form of such a 
regulating mass term M should, ideally, be derived from the dynamics of bubble 
collisions. However it is sufficient for our purposes to know that M ~ for t ~ r c 
and that M — > oo for t ~ r M . This behavior can easily be enforced by the 
"time" -dependence 

M = M Q t a (36) 

where M Q > 0, a > are constants. The effect of the addition of M in to the 
diffusion equation can be clearly seen in the, simple and instructive, case of Ai = M m . 
Then the "massive" heat kernel KM{4>i,4>2',t) gets expressed in terms of the "massless" 
heat kernel i^o(0i,02;^) by 

K M (h, 2 ; *) = Mfa, 2 ; t) e~ tM2D (37) 
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Combining (36) and (37), we observe that the presence of M suppresses super- 
exponentially the contributions of the large t to the heat kernel, which is the de- 
sired result. The basic role of (37) rests in the fact that the "massive" heat kernel 
Km(<I>i, fa', t)M for any generic vacuum manifold M. is expressed as an asymptotic 
expansion of Km(4>i, fa',t) K k in terms of t through [34] 

K M (fa,fa;t) M = K M (fa, fa;t) R k {l + ta 2 (fa,fa) + t 2 a 4 (fa,fa) + . . .} (38) 

Here 02,04,... are polynomials of the Riemann tensor of Ai, its contractions, its 
covariant derivatives, and subsequent combinations of them [37]. We see that, due to 
(38) the super-exponential decay of the large t, contributions still persist for any M. 
as in the case of M, m . The effect of suppressing the large t contribution through M 
can also be seen at the level of the scalar propagator on M. which is [37] 

poo 

Q~ 1 {fa,fa)M = K M (fa,fa;t) M dt (39) 
Jo 

which when combined with (37) gives 

POO 

Q~\fa, fa) M = / K Q {fa, fa; t) M e- tM2D dt (40) 
Jo 

Second, we have to notice that everything we have mentioned so far applies for clas- 
sical or zero-temperature field theories. But it is far more realistic, and thus desirable, 
for this formulation to be applicable when quantum or thermal contributions are taken 
into account. Although these two classes of contributions are quite distinct physically, 
they are handled by very similar methods, within the regime of perturbation theory. 
For this reason from now on, we will refer only to quantum corrections. Quantum 
corrections can be taken into account by replacing the potential energy density V by 
the effective potential V e ff [38] in (1). To illustrate our points we rely, for simplicity, 
on one of the most studied models, the 4 model, whose classical Lagrangian is, for 
a = 1, 2, ... ,k 

Hfa = \(d a fa(d a fa + l -m*tf + A 4 (41) 

Its quantum effective potential V e ff has a loop expansion given by 

v eff = V cl + W 1 + h 2 V 2 + ... (42) 

where V d = ^fa 2 + ^0 4 . The one- loop corrections to V e ff receive contributions 
from the following two mometum integrals, for e > 

* ' ' IfjflV** \ .„ (43) 



4! J (2vr)V -m 2 + ie ' A\ 2 J J (27r) fc (p 2 - m 2 + ie) 
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for the propagator and the four- vertex diagram, respectively [38]. These integrals show 
explicitly how the dimension k of spacetime R fc enters the effective potential. Since 
the vacuum minimizes V e //, the quantum corrections, unlike the classical potential, 
bring an explicit dependence of V e ff on k. This behavior is also in accordance 
with the experimental data: the density of topological defects is observed to depend 
on the spacetime dimension [1]. In the physically important case k — 4 if we carry 
out explicitly these integrals we find that [38] 

mV A0 4 1 2 \ 2 A A0 2 \ A0 2 /3A0 2 2 \\ 

(44) 

We observe that this expression has a divergence for rh — 0. This is due to the 
ultraviolet subtractions in the one loop contribution, which were set up to enforce the 
condition d J^l 1 1^=0 = 0. In case we want to extrapolate the calculation of V e ff to 
rh — 0, we choose a new renormalization scale <p = M and a renormalized coupling 
constant A M at that scale and we obtain the familiar Coleman- Weinberg expression 
[39] 

V - A - 04 + X ^ fin 02 25> l + (45) 
Veff ~ Am 4! + (16^ V M 2 6/ (45) 

We observe that K// has a singularity when In = 2 #. It is very typical for 
V e // to have such singularities. Upon minimization of V e // to determine the vacuum, 
such singularities will persist, will become more numerous and get more complicated 
as the order of the loop expansion increases [38]. As a result, the vacuum M will 
inherit such singularities, so when the quantum corrections are taken into account we 
can no longer assume M. to be a manifold. Although the exact geometric role of 
such singularities is unknown, it is clear that the underlying differentiable structure 
of M. loses its smoothness on them. Since higher orders of perturbation theory keep 
bringing more and more such singularities in the determination of A4, we wonder 
whether it makes any sense to speak about a smooth structure of M at all ! One 
way out of such a difficulty would be to excise, by hand, all the singularities of M. 
This is an iterative and ad hoc process, however, so it is less than satisfying, theoret- 
ically. Another approach, would be to enlarge the structure that M can be allowed 
to have, from that of a manifold to that of a metric space. The advantage of this 
enlargement is that "mild" singularities are already incorporated in the formalism of 
metric spaces [27], [29]. Moreover, there is no a priori requirement of smoothness of 
such a M. In this context, it appears that the metric and measure structures of 
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M., which now, unlike the Riemannian case are disjoint, will keep existing with appro- 
priate modifications, despite the presence of singularities. Therefore, when quantum 
corrections are taken into account we may have to abandon the Riemannian structure 
of M in favor of a metric-measure space structure [40], [42]. In such a case one can 
still define a diffusion equation like the one that gave rise to the geodesic rule, even 
without any smoothness assumptions [41], [42]. We find it very unlikely that in such 
a case the vacuum will maintain its positive sectional curvature (now defined in the 
sense of Alexandrov-Toponogov) [42], as in the classical/non-thermal case. However, 
we expect that relative volume increments (as defined through a Hausdorff measure) 
of M. would have an upper bound, a fact which, in the Riemannian case, corresponds 
to a lower bound on the Ricci curvature [30]. Lower Ricci curvature bounds are neces- 
sary for the stochastic completeness of the diffusion equation in the Riemannian case 
[34] and this may still turn out to be true for the metric-measure spaces of interest [42]. 

Last, we provide a sample of an analytic computation regarding the probability 
of formation of vortices (string-like defects) in a system whose vacuum is RP 2 . An 
example of such a system is furnished by nematic liquid crystals. This calculation can 
also be found in [43]. We view RP 2 as the unit sphere § 2 with opposite points 
identified. Such a vacuum can arise as a result of the symmetry breaking pattern 
5*0(3) — > SO (2) x Z 2 , for instance. To begin with, it is possible to have vortex 
formation in this case, since iti(RP 2 ) = Z 2 . To compute this probability, we have to 
assume that a vortex is formed from the collision of three bubbles, the order parameters 
of which (pi, i — 1,2,3 are independent random variables uniformly distributed, with 
respect to the Lebesgue measure, namely the area, on RP 2 . Without loss of generality, 
we can take 0i located at the north pole of § 2 , which is the double covering of RP 2 . 
The probability we are seeking places the second point 2 at an angular distance of 
at least 7r/2 from 1; since we are interested in a vortex formation. This probability 
is given by the ratio of the area of the spherical cap around the north pole (0i) of 
angular radius 9 over the area of the northern hemisphere. We consider the northern 
hemisphere since a value of <p 2 inside it does not give rise to a vortex, according to 
the geodesic rule and what follows. Such an area ratio equals || = | , for < 9 < it. 
Then <p 3 can be located uniformly anywhere in the hemisphere whose equator is the 
great circle joining 0i with 2 since 3 G RP 2 instead of S 2 . By using spherical 
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coordinates, the probability of formation of a vortex is given by 

-sm9d9 = - (46) 

7T 7T 

This type of integral geometric computation is possible analytically only in very lim- 
ited cases, in which the geometry is so simple as to allow a straightforward enough 
parametrization of M. and subsequent explicit integrations. In all other cases, such 
as almost all computations of the density of higher codimensional defects, a numerical 
estimate seems to be the only feasible way to reach any results. 
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